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Abstract 

In this letter, a joint robust transmit/receive adaptive beamforming for multiple-input multiple- 
output (MIMO) radar based on probability-constrained optimization approach is developed in the case 
of Gaussian and arbitrary distributed mismatch present in both the transmit and receive signal steering 
vectors. A tight lower bound of the probability constraint is also derived by using duality theory. The 
formulated probability-constrained robust beamforming problem is nonconvex and NP-hard. However, 
we reformulate its cost function into a bi-quadratic function while the probability constraint splits 
into transmit and receive parts. Then, a block coordinate descent method based on second-order cone 
programming is developed to address the biconvex problem. Simulation results show an improved 
robustness of the proposed beamforming method as compared to the worst-case and other existing 
state-of-the-art joint transmit/receive robust adaptive beamforming methods for MIMO radar. 
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I. Introduction 

Multiple-input multiple-output (MIMO) radar has beeome the focus of intensive research [[Tl- 
0. MIMO radar with colocated antennas allows to form a virtual array with a larger number 
of virtual antenna elements. It significantly enhances resolution O and parameter identifiability 
[|6l, iQ, allows for direct applicability of adaptive arrays for target detection [[^, dH, [HI, and 
enhances flexibility for transmit beampattern design lfT^ - lfT4]l . 

Robust adaptive beamforming has been widely used to achieve higher resolution capability 
in the traditional frame work of phased receive array [|T5]| - [I^ . The well-known sample ma¬ 
trix inversion (SMI) beamformer llTSll - lIT^ has been recently used in the MIMO radar (joint 
transmit/receive) beamforming context in [|30l, lISTIl . However, the SMI beamformer does not 
provide sufficient robustness against a mismatch between the presumed and actual transmit and 
receive steering vectors. The so-called loaded SMI (LSMI) [|23, [|24l, fl^ and the worst-case 
optimization-based beamformers lfT7l - l(29l can be used, but they may be overly conservative in 
practical applications, because the actual worst operational conditions may occur with a rather 
low probability. Robust adaptive beamforming approaches developed in Il33l - [l36ll for traditional 
phased arrays are then more appropriate since they are designed to guarantee the robustness 
against the signal steering vector mismatch with a certain selected probability. 

In this letter, we propose a joint robust transmit/receive adaptive beamforming for MIMO 
radar using the probability-constrained optimization and provide a mathematical analysis of the 
tight lower bound for the probability constraint by using duality theory for the case of arbitrary 
distributed mismatch. We transfer the primal nonconvex optimization problem into a biconvex 
problem [|371 and address it by using block coordinate descent (BCD) approach [I38l - ll42l . 

II. Signal Model 

Consider a MIMO radar system equipped with transmit and Mr receive antenna elements. 
Both the transmit and receive arrays are assumed to be closely located so that they share the 
same spatial angle of a far-field target. Let 0(f) = [0i(f),... ,0Mt(f)]^ be the waveform vector 
that contains the complex envelopes of orthogonal waveforms which are emitted by different 
transmit antennas, i.e., Jrp(j)i{t)(j)*{t)dt = 6{i — where T is the pulsewidth 

and ^(■) is the Kronecker delta function. Here (■)^ and (■)* stand for the transpose and complex 
conjugate operations, respectively. 
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Assuming that k targets are present, the Mj. x 1 reeeived eomplex observation veetor ean be 
written as 


k 

= Y1 + n(t, r) (1) 

i=l 

where r is the slow time index, i.e., the pulse number, /3j(r) is the refleetion eoeffieient of 
the ith souree with varianee cr^, at(6') and ai.(6*) are the transmit and reeeive steering vectors, 
respectively, and n(f, r) is zero-mean white Gaussian noise. The reflection coefficients l3i{T),i = 
1,... ,k are assumed to remain constant during the whole pulse, but vary independently from 
pulse to pulse. 

By matched filtering the received data to the Mt orthogonal waveforms at the receiving end 
and stacking the individual vector components in one column vector, the MtMj- x 1 virtual data 
vector can be obtained as 


y(r)=vec( / x(f,r)0 {t)dt 

k 


( 2 ) 


= ^ /3i(r)at(6'i) 0 aj.(6'*) + n(r) 

i=\ 

where n(r) is the MtM^ x 1 noise vector whose covariance is given by vec(-) is the 

operator that stacks the columns of a matrix into one column vector, 0 denotes the Kronecker 
product, and (■)^ stands for the Hermitian transpose. 

In contrast to the traditional phased-array radar, the mismatches existing in both the transmit 
and receive steering vectors have to be considered in MIMO radar. Then, the actual transmit 
and receive steering vectors can be modeled as ll^ 


+ Gt, ||et|| ^ St] 


~ A 
Ur = 


— Uj. + Or 




(3) 


where Ot and are unknown complex vectors describing the transmit and receive steering vector 
mismatches, respectively, and H-H denotes the Euclidean norm of a vector. 

In the mismatched case, the actual virtual steering vector can be constructed as 


d = at 0 Uj. = (at -f Ot) 0 (ar + Or) 


= at 0 ar + at 0 Or + et 0 ar + Ot 0 Or. 


(4) 
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The bound on the norm of the virtual steering veetor mismatch can be then found as 

||e|| = ||at 0 Or + et ® ar + Ot ® ej 


^ ||at ® ej + ||at 0 ej + ||at 0 ej (5) 

^ + \/= £• 

It can be seen that the bound @ which is composed of three terms is much larger than that of 
Ot and Or considered separately. Thus, the worst-case approach which uses the loose bound ([^ 
may be too ovely conservative. 


III. Probability-Constrained Optimization 


Since the worst-case approach may be too overly conservative, especially in MIMO radar 
context, we suggest to use the probability-constrained optimization-based approach [|^ - [lT6l . 
The key idea of this approach is to maintain the beamformer distortionless response only 
for operational conditions which occur with a sufficiently high probability rather than for all 
operational conditions corresponding to the uncertainty set. Then the joint transmit/receive robust 
adaptive beamforming problem for MIMO radar can be formulated as 

min w'^Rw 

( 6 ) 

s.t. Pr{w ^(at0ar)^l}^p 

where p is a certain probability value, which can be selected according to the quality of 
service (QoS) requirements, Pr{-} stands for the probability operator, Ot and Or are assumed 

L 

to be random, and R = y(T)y^(r) is the sample covariance matrix used in practical 


r=l 

applications. Here L is the training sample size. 

The general probability-constrained problem (|^ is nonconvex and NP-hard. However, thanks 
to the MIMO radar structure, the joint transmit/receive beamforming vector has the following 
Kronecker form w = u0v. Thus, we can reformulate the cost function of @ into a bi-quadratic 
cost function. 

Using the equalities a 0 b = vec (ba^) and vec (A)'^ • vec (B) = tr we obtain that 

(at 0 dr) = (u^dt) (v^dr). Then the probability constraint in @ can be reformulated as 

Pr{w^ (dt 0 dr) ^ 1} = Pr{|u^dt||v^dr| ^ 1} 


^ Pr{|u^at| ^ 1 n Iv^arl ^ 1} ^ p. 


(7) 
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Since et and ej. are independent identieally distributed (i.i.d.), and using the faet that any funetions 
of independent random variables are statistieally independent ll43l . we ean split the eonstraint 
o into two eonstraints as 

Pr{|u^at| ^ 1 } ^ r/i, Pr{|v'^ar| ^ 1 } ^ 772- ( 8 ) 


The probability values 7]i and r ]2 are ehosen separately at the transmitter and reeeiver so that 
^ 1^2 = P and ^ 1- Note that rji and t ]2 are preseleeted based on the robustness requirements 
to mismatehes at the transmitter and reeeiver. Thus, rji and r ]2 are not the optimization variables. 
For example, pi = P 2 = y/p satisfies the eonditions. 

Let us obtain a simplified approximate form of the probability eonstraints @. If |u^at| > 
|u^et| and |v'^ai.| > |v^er|, i.e., if the steering veetor mismatehes are reasonably small, then 


from the triangle inequality it follows that 


|v^ (a^ + Or 


lu'^at 

— 

|u^et| 

|vV| 


v^e,| 


(9) 


Using @ and @, the problem (|^ ean be reformulated as 

\H 


mm 

U.V 


u 


R u 


s.t. Pr {lu'^Otl ^ |u^at| — l} ^ Pi 
Pr jlv^erl ^ |v'^ar| — l} ^ p 2 - 

The problem ( fT^ beeomes mathematieally traetable if we additionally assume a speeifie 
analytie form for the probability operator Pr{-} and make some approximations. In the sequel, 
we will eonsider two praetieally important oases eorresponding to two different assumptions on 
the probability density funotion (pdf) of the transmit and reoeive steering veetor mismatehes Ot 
and Of. 


A. Mismatches With Gaussian Distribution 

Consider the ease of a oomplex zero-mean symmetrio Gaussian distribution for and Or, i.e., 

et~CAr(OM„Ct), er~CA^(OM„Cr) (11) 

where Om* and denote the Mt x 1 and Mj. x 1 veetor of zeros, respeetively, Ct and Cr 
eapture the seeond-order statistios of the uneertainties in the transmit and reeeive steering veetors, 
respeetively. . 
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We consider the receive steering vector mismatch and note that the same considerations apply 
to the transmit steering vector mismatch. It is easy to show that the random variable v^Or has the 
complex Gaussian distribution, that is, v^Or ~ CA/^ (Om^, v^Crv), and its real and imaginary 
parts are real i.i.d. Gaussian. 

Let us use the fact that if x and y are two real i.i.d. zero mean Gaussian random variables 
with the variance cr^, then z = is Rayleigh-distributed with the cumulative density 

function (cdf) given as F{z) = 1 — Using this fact, the probability constraint @ at the 

receiver can be written as 


Pr {Iv^erl ^ |v^ar| — l} 


= 1 — exp 


(|v^ar| - 1)^ 


( 12 ) 




m- 


The inequality in the second line of ( fT^ can be equivalently rewritten as 

cIv 


v^arl — \ lln 


1 -??2 


^ 1 . 


(13) 


Observing that the cost function in ( fTO] ) is unchanged when u and v undergo an arbitrary 
phase rotation, the problem (fTO]) can be further rewritten as 


min (u 0 v)^R(u®v) 

U,V 


s.t. u^at - 7ii 

C|u 

^1, 

Xm(u^at) 

= 0 

(14) 

v^aj. - 712 

cIv 

^ 1, 

Xm(v^ai.) 

= 0 



where yn = ^/n((l and 712 = ^/n((l - 72 )"^)- 

The problem ( [T?] ) has separate constraints for u and v, but coupled objective. For such type of 
problems of minimizing a continuous function of several blocks of variables, the block coordinate 
descent (BCD) methods are widely used [|^ - Il4^ . At each iteration of BCD, a single block of 
variables is optimized, while the remaining variables are held fixed. Thus, we develop here a 
BCD type method for addressing the problem ( [T4| ). It can be seen that if one of the vectors u or 
V is fixed, the cost function of the problem ( [T4| ) can be transformed into a quadratic function with 
respect to the other vector. Hence, it can be solved based on the second-order cone programming 
(SOCP) using, for example, the CVX toolbox Il44ll . 
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Algorithm 1 : Interative SOCP Method 
1 : Given: initial value uq = /c = 0, toleranee S 

Repeat: k = k + 1 

2 : Solve ([T5]), obtain := 

3: Solve ( [T^ , obtain u^. := 


Until 


4: 


opt I 

H-M 


Return 


< 6 and 


opt I 

11 opt I I 

r^o-i 


< 5 


5: Output: u°P* and 


Let us fix u and ehoose It is obvious that there is a sealing determinaey between 

u and V. Then, ([T?]) boils down to 


min (ufc (g) v)^ R (ufc 0 v) 


s.t. (^uf at - 7ii 

Xm(v^ai.) = 0 




v^ar - 7 i2 


civ 


> 1 


(15) 


By solving ( [T5] ), we obtain the optimal solution Then we fix v and ehoose Vfc = 

The problem ([T?]) boils down respeetively to 


min (u (g) Vfc)'^ R (u (g) Vfc) 


U 

s.t. 



at - 7ii 




ar 


— 7l2 


cl 


Vfc 


^ 1 


(16) 


Xm(u^at) = 0. 

The overall algorithm for addressing ( fT?] ) is then given as in Algorithm 1. The eonvergenee of 
a large elass of BCD methods is investigated in ll4Tll . 


B. Mismatch With Arbitrary Distribution 

Consider now the ease when the transmit and reeeive steering veetor mismatehes are arbitrary 
distributed and only the first and seeond-order statisties are known. 

Using the Chebyshev inequality whieh states that for any zero-mean random variable r with 

2 

varianee cr^ and positive real number a, we obtain Pr{|r| ^ a} ^ Then the probability 
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constraint @ at the receiver can be expressed as 


Pr{|vV|^|vV|-l} = 1 - 


(|v^ar| - 1)2 

Further, ( fTT] ) can be rewritten at both the transmit and receive sides as 


^ r} 2 - 


|u^at| -721 




^ 1, |v^ar|-722 

-1 


Civ 


^ 1 


(17) 


(18) 


where 721 = (-^/l — 71 ) , and 722 = (\/l — 72 ) • We can see that the difference between the 

cases of the mismatch with arbitrary and Gaussian distribution is in the coefficients 711,712 and 
721,722 only. Thus, the problem in the case of the mismatch with arbitrary distribution can be 
addressed in the same way as before. 

IV. Tight Lower Bound for The Probability Constraint Using Duality Theory 

For simplicity, let us eliminate the subscripts in our notations and instead use simply the 
notations a, e, and C. The tight lower bound problem is 


min Pr{|v^ (a + e)| ^ 1} 
/(e) 

where /(e) is the pdf of e. It can be reformulated as 

min Pr{ 5 f(e) ^1} 

/(e) 


(19) 


( 20 ) 


s.t. E{ee^} = C 

/ f{e)de = 1, E{e} = 0 
J e 

where g{e) = |v^(a + e)|2. To find the dual problem for pO] ), we first introduce the Lagrangian 
function as 

L(/(e), /i, T}, T) = Pr{^(e) ^ 1} + /r(l - f f{e)de) 

J e 

+ r/^(0 - E{e}) + tr(T^(C - E{ee^})) 

= /i + tr(T^C) 

+ [ [1 - (yu + + tr(T^ee'^))]/(e)de 




( 21 ) 


[0 - (/r + ri^e + tr(T"ee'^))]/(e)(ii 


H _ H\ 


G(e)s£l 
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where /i, rj, and T are the Lagrange multipliers and T = T^. With the implieit pdf eonstraint 


/(e) ^ 0, the Lagrange dual funetion of the problem (20) is given as 


= min L(/(e),/i, r/, T). 

/(e)^0 


( 22 ) 


Note that the minimum of the first integral in ( |2T] ) with the nonnegative pdf /(e) eonstrain is 
zero if /i + rj^e + tr(T^ee^) ^ 1, otherwise the minimum is unbounded below. Similarly, the 
minimum of the seeond integral with the nonnegative pdf /(e) eonstrain is zero if /i + rj^e + 
tr(T'^ee^) ^ 0, otherwise the minimum is unbounded below. It is easy to see that the eondition 
/x + ? 7 ^e + tr(T^ee'^) ^ 0 must hold for both g(e) ^ 1 and g{e) ^ 1, i.e., \fg{e), and therefore, 
Ve, while the eondition jj, + rj^e + tr(T'^ee^) ^ 0 must hold for yg{e) ^ 1. Then the dual 


problem to (20) ean be formulated as 


max /i + tr(T"C) 

s.t. fji + rj^e + tr(T^ee'^) ^ 1, Ve 

fji + r/^e + tr(T^ee'^) ^ 0, ^gie) ^ 1 

rp _ rj^H 


(23) 


Let us define 


~ A 

e = 


,C^ 


C 0 

0 ^ 1 


z ^ 










vv 


H 


vv'^a 


a^vv'^ a-^vv^a — 1 


Thus, the seeond eonstrain in (23) ean be rewritten as e^Ze ^ 0, Ve^ Ae ^ 0. Further, aeeording 


to S-Lemma it beeomes Z — AA ^ 0, VA ^ 0. Then the problem ( [23| ) in eompaet form ean be 
expressed as 


max tr(ZC) 


s.t. Z — 


0 0 

0^ 1 


V 0, Z — z 


H 


(24) 


Z - AA ^ 0, VA ^ 0. 
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Since the Lagrange dual problem yields the lower bound of the primal problem, we ean finally 
state that 

min Pr{ 5 f(e) ^ 1 } = max tr(ZC). (25) 

/(e) Z 


V. Simulation Results 


We assume a uniform linear array (ULA) of 10 antenna elements that is used for both 
transmitting and reeeiving, i.e., Mt = 10 and Mr = 10. Antenna elements are spaeed half a 
wavelength apart from eaeh other. The plane-wave target impinges on the array from Og = 3° 
and two interfering sourees arrive from direetions On = 30° and 6*12 = 50°, respeetively. The 
interferenee-to-noise ratio (INR) is assumed to be 20 dB for both interferenees, and 100 Monte- 
Carlo runs are used to obtain eaeh point in our simulations eurves. 

Consider the seenario with Rieean propagation medium where the mismateh veetors Ot and 
Or are modelled as 

N 


e = 


a 

^MN 


J2e^^M0o + 0n) 

n=l 


where M stands for Mt or M^, a is the power of seattered nonline-of-sight (NLOS) eomponents, 
N is the number of NLOS eomponents, ^jJn is the phase shift parameter of the nth NLOS 
eomponent, and 6n is angular shift. The parameters 9n and ipn, n = 1,... ,N are independently 
and uniformly drawn in eaeh simulation run from [—2.5°,2.5°] and [0,27r), respeetively. The 
values = 10 and cx^ = 0.3M are taken. 

For eomparison, e = 9 is ehosen for the worst-ease robust beamformer [(32l. For the LSMI 
beamformer, the fixed diagonal loading parameter 7 = 10 is ehosen. For the proposed probability- 
eonstrained optimization based joint transmit/reeeive robust adaptive beamformer, p = 0.9 while 
rji = 0.93 and 72 = p/Vi = 0.9677 are taken. 

The output SINRs versus SNR are shown in Fig. We ean see that the proposed method 
based on the probability-eonstrained optimization for the ease of arbitrary distributed mismatehes 
has the best performanee among all the teehniques tested. These improvement is espeeially 
remarkable at high SNRs. However, the performanee of the probability-eonstrained optimization 
based method with Gaussian mismateh distribution is worse than that for the worst-ease robust 
beamformer. It ean be explained by the faet that the aetual mismateh eorresponding to the 
eonsidered Rieean seenario is not Gaussian. 
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Fig. 1. Output SINK versus SNR for K=100 and INR=20dB. 


VI. Conclusion 

A joint transmit/receive robust adaptive beamforming method for MIMO radar has been 
developed based on the probability-eonstrained optimization. Specifically, we have considered 
the mismatches in the desired signal steering vectors at both transmit and receive arrays to 
be random. The original probability-constrained optimization problem has been converted into 
a biconvex problem and addressed by using BCD approach. The proposed technique offers a 
better performance than several state-of-the-art counter parts. 
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